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SELF- CONSISTENT  NONLINEAR  SLOW- TIME  SCALE  FORMULATION  AND  SIMULATION  OF 
OVERMODED  GYROTRON  OSCILLATORS  AND  AMPLIFIERS 

1.  INTRODUCTION 

Hie  gyrotron  (or,  electron  cyclotron  maser)  instability  arises  from 
the  interaction  of  a  TE  waveguide  (or  cavity)  mode  with  a  relativistic 
electron  beam,  which  drifts  along  and  rotates  about  an  externally  applied 
static  magnetic  field.  This  instability  occurs  for  sufficiently  close 
synchronization  of  the  relativistic  electron  gyrofrequency  in  the  external 
field  to  the  wave  frequency,  Doppler  shifted  in  the  electron  frame,  since 
it  is  then  that  the  wave  electric  field  acts  on  the  beam  electrons  co¬ 
herently,  causing  changes  in  their  relativistic  factor,  y,  depending  on 
the  electron  phase  relative  to  the  wave,  resulting  in  phase  slippage  and 
bunching,  with  the  corresponding  current  driving  in  turn  the  wave.  Volu¬ 
minous  theoretical  and  experimental  investigations  of  this  instability 
(1-12  and  references  therein)  have  established  the  gyrotron  as  a  very 
potent  source  for  microwave  generation,  with  inherent  efficiencies 
approaching  50%.  Further  improvements  to  this  value  for  the  efficiency 
are  possible  by  introducing  axial  gradients  either  in  the  diameter  of  the 
oscillator  (13,  14),  hence  the  axial  wavenumber,  or  in  the  guiding  field 
(15-17). 

In  view  of  the  impressive  performance  of  the  gyrotron  in  converting 
the  kinetic  energy  of  the  beam  electrons  to  electromagnetic  energy,  it  is 
the  prime  candidate  to  be  used  as  the  radiation  source  for  the  heating  of 
tokamak  plasmas  at  the  electron  cyclotron  frequency.  For  the  parameters 
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currently  considered,  e.g.,  for  a  confining  field  of  the  order  of  50  kG, 
requiring  a  radiation  wavelength  of  the  order  of  2  ran,  the  obvious  size 
and  power  limitations  require  an  operation  of  the  gyrotron  not  at  the 
fundamental  harmonic.  Accordingly,  it  is  required  that  the  operation  of 
an  overmoded  gyrotron  cavity  be  investigated,  to  assess  the  performance 
of  the  desired  mode  and  the  potential  of  destructive  interference  by  the 
competing  modes. 

The  problem  of  multimode  interactions  in  an  electron  cyclotron  maser 
has  been  addressed  analytically  in  the  Soviet  literature  (IP- 20).  In  their 
analysis,  the  coupled  system  of  equations  for  the  eigenmodes  and  the  beam 
distribution  is  presented,  and  the  coupling  coefficients  are  expressed  in 
terms  of  series  involving  the  even  powers  of  the  field  amplitudes.  From 
the  consideration  of  the  coefficients  of  these  expansions,  distinction 
between  soft  and  hard  coupling  is  made  and  general  conditions  are  given 
to  distinguish  between  stable  and  unstable  equilibria.  This  approach  is 
mathematically  involved  to  an  extent  that  obscures  the  physics  of  the  inter¬ 
action,  and  its  applicability  appears  limited  to  situations  involving 
restrictions,  such  as,  a  low  amplitude  for  all  but  one  mode,  or  small 
beat  frequencies  compared  to  the  inverse  of  the  electron  transit  time 
through  the  resonator,  or  a  premodulated  electron  beam,  or  an  identical 
axial  dependence  of  the  modes. 

An  alternate  approach  has  been  recently  initiated  at  NRL  (21),  in 
which  the  motion  of  test  particles  is  followed  in  cavity  waves  with  constant 
amplitude  and  frequency,  the  latter  taken  equal  to  the  empty  cavity  eigen 
frequency.  Equating  the  energy  lost  by  the  particles  to  each  mode  to  the 
diffraction  losses  gives  the  values  of  the  quality  factors,  Q,  consistent 


with  equilibrium,  and,  by  iteration,  the  equilibrium  amplitudes  (and 
efficiencies)  for  given  values  of  Q  (and  of  the  beam  power) .  A  presentation 
of  these  results  in  the  form  of  contour  plots  appears  to  be  very  useful 
for  the  determination  of  the  operating  conditions.  However,  these  results 
are  liable  to  reservations  regarding  the  assumption  of  constant  (and  not 
oscillating  at  the  beat  frequency)  amplitudes  and  frequencies,  and  that 
the  latter  is  not  shifted  by  the  action  of  the  beam. 

The  purpose  of  this  work  has  been  to  develop  and  apply  a  self- 
consistent  formulation  for  the  multi-mode  problem,  which  is  simultaneously 
adoptable  for  efficient  numerical  simulations  and  applicable  for  theoretical 
considerations.  After  a  brief  description  of  the  physical  model  in  Sec.  2, 
the  theoretical  method  is  discussed  in  Sec.  3,  while  the  resulting  equations 
are  presented  in  Sec.  4,  along  with  a  discussion  on  the  algebraic  manip¬ 
ulations  that  led  to  their  derivation.  In  order  to  establish  the  correctness 
of  the  equations  and  to  demonstrate  their  applicability  for  an  analytical 
treatment,  they  are  used  in  Sec.  5  to  rederive  and  generalize  the  equations 
pertinent  to  the  evolution  in  a  waveguide,  and  this  evolution  is  shown  in 
Sec.  6  to  be  accurately  described  by  the  algorithm  based  on  these  equations. 
Some  consideration  of  introductory  nature  is  given  in  Sec.  7,  on  the  inter¬ 
action  responsible  for  the  efficiency  enhancement  in  a  nonuniform  external 
field,  while  examples  on  the  application  to  an  overmoded  cavity  are 
presented  in  Sec.  8.  Finally  Sec.  9  summarizes  our  results,  compares  our 
method  to  the  test  particle  approach,  and  indicates  various  possible 
extensions  to  this  work. 
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2.  DEFINITIONS 

This  section  serves  the  purpose  of  defining  the  physical  model  we  are 
considering  and  the  nomenclature  for  the  quantities  of  interest. 

We  consider  a  rectangular  cavity  with  perfectly  conducting  walls.  The 
transverse  dimensions  of  the  cavity  are  xmax  and  ymax»  and  the  length  is  z(nax 
(Fig.  1).  The  choice  of  the  rectangular  geometry  was  for  the  purpose  of 
simplifying  the  algebra,  while  the  physics  of  the  interaction  is  not  obscured. 
The  results  obtained  with  this  geometry  can  be  directly  applied  to  a  coaxial 
cavity  with  large  aspect  ratio,  while  for  a  general  cylindrical  geometry  an 
appropriate  conversion  of  either  the  results  or  the  formalism  is  necessary. 

An  external  magnetic  field  is  present  inside  the  cavity.  The  main 

component  of  this  field  is  B^ex^,  along  the  axis  of  the  cavity,  however,  a 

weak  nonuniformity  is  allowed,  i.e. ,  9B^ex)/3z  and  the  transverse  components, 

B^ex^  and  B^ex\  are  small  but  not  necessarily  zero.  The  introduction  of  this 
x  y 

nonuniformity  is  for  the  purpose  of  assessing  the  associated  efficiency  enhance 

ment  (15,  16),  similar  to  the  corresponding  behavior  in  the  case  of  the 

Whistler  mode  (22,  23).  For  normalization  purposes,  it  is  appropriate  to 

express  the  external  field  components  by  their  nonrelati vistic  electron 

gyrofrequencies,  fix,  siw  and  nz,  defined  by  Slz  =  eBzexVmec,  with  similar 

definitions  for  n  and  SI.,  where  -  e  and  m  are  the  electron  charge  and  rest 
x  y  c 

mass,  respectively,  while  c  is  the  speed  of  light. 

An  electron  beam  is  continuously  injected  into  the  cavity  at  z  =  0 
and  propagates  self-consistently  with  the  external  magnetic  field  and  the 
cavity  fields,  defined  below.  The  distribution  function  of  the  electrons  is 
defined  as 
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(1) 


^max^ 

f(r,  v,  t)  =  L<S(v  -  j^(t))  6(r  -.^.(t)). 


m 


A 


where  £.(t)  and  v,.(t)  are  the  instantaneous  position  and  velocity  vectors  of 

the  i-th  particle,  obtained  from  its  equations  of  motion  in  terms  of  its 

entry  properties,  and  imax(t)  is  the  number  of  particles  inside  the  cavity 

at  time  t.  No  assumption  is  made  with  regard  to  the  form  of  the  distribution 

function  at  z  =  0,  hence  it  is  possible  to  consider  the  effect  associated  with 

a  position  or  velocity  spread  of  the  electrons  or  with  any  slow  time  dependence 

of  their  bulk  entry  properties  (e.g.,  velocity  or  density).  It  is  convenient 

to  introduce  the  instantaneous  plasma  frequency,  Wp(t),  of  the  beam,  defined  by 

<*£  =  (4ire2i‘  x  w  y„av  zma„).  It  is  assumed  that  the  beam  density  and 

P  max'  '  e  max  •'max  max'  J 

axial  velocity  inside  the  cavity  are  not  too  much  different  from  those  outside, 
so  that  any  effects  on  the  beam  distribution  by  its  self-fields  can  be 
adequately  accounted  for  by  an  appropriate  expression  of  the  distribution 
function  at  the  entry  point.  However,  the  exact  (i.e.,  self-consistent) 
particle  density  and  axial  velocities  inside  the  cavity  will  be  considered  for 
the  computation  of  the  beam  effect  on  the  cavity  fields.  The  transverse 
components  of  the  position,  (x,y),  and  velocity  (v¥,  v  ),  of  the  electrons  are 

a  y 

expressed  by  the  position  of  the  guiding  center,  (x,  y),  and  the  amplitude  and 
phase  of  the  velocity,  (v^,  ip),  defined  by 


x  =  x  +  p  sin  ip,  y  =  y  -  p  cos  ip  , 

(2) 

vx  =  VA  cos  ip  ,  vy  =  Vx  sin  Ip  , 

as  shown  in  Fig.  2,  where  p  =  yv  /n  is  the  Larmor  radius  and 

it 

y  =  c/(c2  -  v2  -  v2)*5  is  the  relativistic  factor.  Obviously,  this  transfor¬ 
mation  is  suggested  by  the  fact  that  the  principal  motion  of  the  electrons 
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/  Av\ 

is  that  under  the  action  of  the  axial  external  field  component,  B^  '. 

The  fields  of  the  cavity  are  expanded  in  terms  of  the  vacuum  eigenmodes. 
Since  the  gyrotron  instability  involves  the  TE  modes,  only  these  modes  are 
included  in  the  expansion.  The  (n,  m,  £)  mode  is  characterized  by  the  con¬ 
stant  wavenumbers  ( kn ,  km,  k£)  corresponding  to  the  three  Cartesian  axes  and 
defined  by  kn  =  nir/x^  (n  =  0,1,2,...),  km  =  rmr/ymax  (m  =  0,1,2...)  with 
the  restriction  n  +  m  f  0,  and  k^  =  i/rr/zmax(i,  =  1,2,3...).  The  transverse 

wavenumber  is  defined  by  k  _  =  (k*  +  k*)*5  and  the  total  wavenumber  is 

nm  v  n  m 

k _ _  =  (k*  +  k*)^.  The  4  me  evolution  of  each  mode  is  described  by  the  varia- 

nmx  nm  %' 

tion  of  the  amplitude,  BnniA(t),  or  its  electron  gyrof requency ,  nnmA=  eBnm£/mc’ 
and  the  phase,  4>nm£(t),  or  the  frequency,  <i>nin£  =  d<J>nm£/dt*  In  addition,  we  use 

the  growth  rate,  da^/dt,  and  the  frequency  shift,  Anmt  =  ^  - 

knm£,  c’  as  wel1  as  the  factor,  QnmJl»  which  is  assumed  to  be  much  larger 

than  unity.  Using  the  above  definitions  and  the  shorthand  notations  Cn  =  cos 
(knx) >  sn  =  sin(knx),  Cm  =  cosfk^),  Sm  =  sinfk^),  C£  =  cosfk^z)  and 
=  sinfk^z),  the  cavity  fields  are  given  by 


B^w)  =  £  cn  C  S.  B  #  sin  d>  0  , 
z  n  m  x,  nmx  nmx 

nmx 


B(w)  =  .  y  knlct 


S„  Cm  C0  Bnm0  sin  4.nmjl  , 


k*  n  m  i  nm£ 
nmx  nm 


B 


k  k 

,  = '  S,  €r  c"  s"  c*  s1n  ■ 

nm£  nm 


(3) 


E^W)  =  0  , 


:(W) 


m 


Bm  B0  lf+  (Bnmff  S^n  $  nm£^  * 


‘X  ZHo  cknm  n  *  3t  v  nm£ 

nmx  nm 
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m 


Ey^  =  "  ?»  ck*"  Sn  Cm  si  3t  (Bnm Jt  s1n  <t,nmA^  ' 

J  nmx  nm 

for  the  components  of  the  magnetic  induction  and  electric  field  along  the  z, 
x  and  y  axes,  respectively.  The  summation  formally  extends  over  the  triple 
infinity  of  eigenmodes  supported  by  the  cavity,  however,  when  our  resultant 
equations  are  to  be  used,  it  suffices  to  sum  only  over  those  modes  whose  forward 
or  backward  component  is  adequately  resonant  with  the  beam. 

The  quantity  of  interest  in  oscillators  is  the  efficiency,  nnmA»  of  con¬ 
version  of  beam  power  to  radiation  power  of  the  mode  (n,  m,Jl).  Calculating 
the  field  energy  from  Eq.  (3),  multiplying  by  Wnm£/Qnm£  and  dividing  the 
result,  the  radiation  power,  by  the  beam  power  yields  the  formula  for  the 
effi  ci  ency, 

-  1  +  ^nnuo  ^nmt  zmax  /Mnmt\  /nnmA  ,  (4) 

n"“_  16  vk™‘/vw 

where  =  1 ,  if  either  n  =  0  or  m  =  0,  and  =  0,  otherwise.  The 

nm,o  nm,o 

efficiency  in  Eq.  (4)  is  calculated  in  terms  of  the  instantaneous  values  of 
“nmJ,  and  finm$,  and  entry  va^ues  Y0»  vZQ  and  wpo  of  the  beam,  which  may  be 
slowly  time  dependent,  if  a  beam  with  variable  properties  is  employed. 
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3.  THEORETICAL  METHOD 

In  principle,  the  solution  to  the  problem  involves  the  solution  of  the 
wave  equation  for  the  cavity  fields,  e.g. , 

(5!  -  Jy  Hr)  B(w)  -  -  ^  (V  x  0)  ,  (5a) 

with  the  current  density,  J^r.t)  =  -  bJ d v  v  f(r,  v,  t),  obtained  from  the 
solution  of  the  equations 

IfX’s  • 

(5b) 

[ce(“)*»x(8<»)+b(->,]  , 

for  the  relativistic  motion  of  the  beam  electron  under  the  joint  action  of  the 
cavity  and  external  fields.  The  above  equations,  together  with  Faraday's 
law,  c  V  x  E^  +  3B^w^/3t,  form  a  complete  system,  whose  solution  can  be 
Fourier  analyzed  in  the  eigenmodes,  Eq.  (3),  to  give  the  evolution  of  the 
spectrum.  Such  an  approach,  although  straightforward  conceptually,  has 
limitations  associated  with  its  numerical  implementation,  since  it  involves 
a  large  number  of  arithmetic  operations.  In  addition,  for  a  theoretical  in¬ 
vestigation  of  the  physics  of  the  multi -mode  interaction,  this  approach  is 
not  useful,  since  its  equations  are  not  formulated  in  the  time  scale  of  the 
energy  transfer  from  the  beam  electrons  to  the  fields. 

The  aim  of  this  work  is  to  develop  a  formulation  which  Is  both  efficient 
for  numerical  applications  and  useful  for  theoretical  considerations,  without 
sacrificing  the  rigor  of  self-consistency.  The  central  assumption  in  our 
approach  is  that  there  exists  a  time  Interval,  At,  during  which  the  energy 

10 


transfer  between  the  electrons  and  the  fields,  and  between  the  parallel  and 
perpendicular  components  of  the  electron  motion,  is  small,  while  At  is  not 
necessarily  small  compared  to  the  gyroperiod  or  the  wave  period.  For  a 
numerical  application  of  our  equations.  At  will  be  the  time  step.  Then, 
according  to  our  assumption,  the  quantities  wnmjl»  ^nmA»  *-j,  y.,-.  and  vz^ 
will  be  constant  to  zero  order  during  the  time  interval  At,  while  the 
quantities  ^  and  ^  will  vary  linearly  with  time, 

^nm£^  “  ^nm£^o^  +  ^  “Vimt^o^  * 

Zi(t)  =  Zi(t0)  +  (t  -  tQ)  vz.(t0)  ,  (6) 

^(t)  =  ^(t0)  +  (t  -  tQ)  flz(2i(t0))/Yi(t0)  , 


where  tQ  is  the  time  at  the  beginning  of  the  time  step  and  t  -  tQ  $  At.  These 
zero  order  solutions  can  be  introduced  into  the  differential  equations  to 
obtain  the  first  order  equations  for  the  electron  propagation  and  wave  evolu¬ 
tion  at  the  time  scale  of  energy  transfer. 

Formally,  the  above  assumption  can  be  described  quantitatively  by  the 
strong  inequality 


kc j  >>  Tj  A*  n 


V  ^  3z 


where  the  subscript  triad  nm£  has  been  suppressed  and  e+  =  0+mnA  j  - 

“nmJl  +  k^  vzi  is  the  rate  of  change  of  the  particle  phase,  i|»^,  relative  to 
the  wave  phase  in  the  particle  rest  frame,  <J>nmA  +  k^.  Using  this  inequality. 


’~*T< 

•V 


each  quantity  of  interest,  say  q,  where  q  stands  for  x,  y,  z,  v^,  vz,  ^ •  r 
and  A  or  any  combination  of  them  (like  y)»  can  be  associated  with  a 
differential  equation,  which  is  expanded  in  the  form 


£  A..  sin(a,t  +  pj  , 

j=o  J  J 


(8) 


where  the  amplitudes  are  required  to  satisfy  A^j  »  A^+1  j  and  A..J  -  A^  j+1, 
and  the  frequencies  o^.  «  Oj+p  with  p^  being  an  appropriate  phase.  The 
quantities  A^,  and  p^  are  assumed  to  vary  slowly  during  the  energy  trans¬ 
fer  time,  i.e. ,  to  be  constant  during  At.  The  zero  order  term,  AQ0,  represents 
the  unperturbed  cyclotron  motion  of  the  electrons  in  a  uniform  external  field 
and  the  steady  state  propagation  of  the  cavity  modes,  while  the  terms  under  the 
summation  represent  the  effects  to  first  and  higher  order  (given  by  i)  of  the 
nonuniformity  and  of  wave-particle  coupling,  and  the  slow  (j=0)  and  fast  (to 
all  degrees,  j)  phase  variations.  The  slowest  phase  variation  corresponds  to 
guiding  center  motion  in  the  nonuniform  external  field  (aQ  =  0)  and  to  resonant 

interaction,  or  =  e.. 

o  + 


Clearly,  there  is  no  need  to  evaluate  and  use  many  terms  in  the  expansion 
of  Eq.  (8).  A  truncation. 


Bt  =  Aoo  +  A10  sin(V  +  Po}  »  (9> 

is  expected  to  be  sufficient,  since  the  contribution  of  the  higher  order  terms, 
i  >  1,  is  negligible,  (A^/A00)  ~  (ak/Aoo^  «  1,  and  the  effects  of  the  fast 
terms,  j  >  0,  phase  mixes  to  zero  during  the  energy  transfer  time,  in  view  of 
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cu/r  ~  ju/r  »  1.  In  this  sense,  our  approach  is  conceptually  identical 
to  that  already  employed  for  electrostatic  and  whistler  waves  (24) . 

A  final  remark  is  associated  with  the  choice  of  time,  t,  rather  than 
the  axial  position,  z,  as  independent  variable.  Certainly,  the  conversion 
of  the  method  and  the  equations  is  straightforward,  by  the  simple  trans¬ 
formation  dt  =  dz/v .  The  reason  behind  this  choice  has  been  to  guard 

it 

against  the  possibility  that  the  parameters  of  some  electrons  become  multi¬ 
valued  functions  of  z,  which  could  well  be  the  case  especially  due  to  the 
nonuniformity,  and  to  avoid  complications  in  the  evaluation  of  the  wave 
properties ,  which  are  necessarily  time  dependent  quantities  and  would  be 
obtainable  by  an  appropriate  interpolation  in  the  electron  dynamics. 


4.  SLOW-TIME-SCALE  EQUATIONS 

In  this  section  we  present  the  slow-time-scale  equations  for  the 
propagation  of  the  electrons  and  the  evolution  of  the  eigenmodes.  Our  approach 
has  already  been  discussed  in  the  previous  section.  Since  the  procedure 
involves  a  lengthy  series  of  calculations  which,  however,  are  of  elementary 
nature,  we  will  limit  ourselves  to  the  presentation  of  the  final  results,  with 
only  a  discussion  of  the  steps  used  in  obtaining  them.  (See  also  Ref.  25.) 

The  slow-time-scale  equations  for  the  electron  dynamics  are: 

li  ■  •  \  v  *r 


-  i  E  Lj  £  ^kmVm  slnc*  +  Wm  cosc*'l  • 

c  nms,  L  nm  +  J 


H  1 

3t  <*vx>  ■  2  V  3T 


•  *  «  h  ^  1  (’  '  !r)k"6"1"  SinCi  +  k"5"E”  “Hi’] 

•ft1*  **  Lt’  JStt?S  ?(k"£"5"'  ""'I  +  k"S"Cm  C0SCi)]  1 


Ht  *  '  7  (S:  + 


?  2  [o,3p-Z(±  i)  (i  -  TpT/Bz) 

c  nmi.  l  YKnm  +  V  Knm  p  / 

x  ( kmC„S_  cose.  -  k„ SC_  sine.)  I  , 
’  m  n  m  +_  n  n  m  —  I 


(10) 
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In  the  above  equations, a  superscripted  bar  denotes  evaluation  at  the  guiding 
center,  i.e.,  Cn  =  cos(knx)  etc.,  and  Jj  are  the  Bessel  function  of  first 
kind  and  first  order  and  its  derivative,  both  evaluated  at  knn)p,  where 
p  =  yv  /fiz  is  the  Larmor  radius,  and  e+  =  V  -  4>nmJl  +  k^z  is  the  electron  phase 
relative  to  the  forward/backward  wave  and  0+  =  dc+/dt  is  the  corresponding 
relative  frequency.  For  numerical  purposes,  these  equations  can  be  Integrated 
analytically  over  the  time  step,  At,  since  during  that  time  all  quantities  of 
the  right-hand  side  can  be  considered  to  be  constant,  except  for  c+,  for  which 


a  variation  dc+/dt  =  C+(tQ)  +  (t  -  tQ)  0+  (tQ)  is  sufficient,  in  terms  of  the 
corresponding  values  at  tQ,  the  beginning  of  the  time  step.  These  integrated 
expressions  can  be  used  to  evaluate  the  integrals  of  dy/dt  and  d(yvz)/dt,  which 
appear  on  the  right  hand-side  of  Eq.  (10),  thus  resulting  in  an  algebraic 
single-step  algorithm.  Although  this  algorithm  includes  linear  equations  in 
a  uniform  external  field,  a  successive  application  of  these  equations  will 
give  the  nonlinear  nonuniform  final  solution,  since  the  nonlinearities  and 
nonuniformities  are  accounted  for  by  the  reinitialization  after  each  time  step. 

Let  us  present  now  a  few  comments  on  the  derivation  of  Eq.  (10).  We 
have  started  with  the  equation  of  motion,  Eq.  (5b),  in  which  the  external  and 
wave  fields  are  treated  separately,  the  latter  expressed  by  Eq.  (3)  and  the 
transverse  position  and  velocity  are  expressed  by  Eq.  (2),  resulting  ‘n  the 
evaluation  of  certain  quantities  of  (x,y),  as  indicated  in  Eq.  (10)  by  the 
superscripted  bars.  With  regard  to  the  contribution  of  the  cavity  fields, 
expressions  similar  to  cos(knmp  cosi^)  are  expanded  in  Fourier 
series,  with  the  coefficients  obtained  from  the  integral  representation  of  the 
Bessel  functions, 

it 

W1  *  1  /  “5<kn»p  s1"*  '  n'rtd*  ’  tn) 

0 

as  well  as  their  symmetry  relations.  The  products  of  trigonometric  functions 
are  converted  to  sums  whose  argument  changes  at  the  frequency  snz/y  +  u>  +  k^vz. 
We  consider  only  interactions  at  the  fundamental  (s=l)  and  elimination  of  the 
fast  terms,  nz/y  +  u  +  kjVz,  leaves  only  the  two  slow  relative  frequencies, 
e+,  and  the  corresponding  phase  angles,  c+,  in  terms  of  which  Eq.  (10)  Is 
expressed.  A  similar  application  of  Eqs.  (5)  and  (2)  on  the  external  field 
contribution,  in  which  terms  varying  at  the  frequencies  nz/y  and  2«z/y  are 
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omitted,  gives  the  expected  acceleration  due  to  the  nonuni formity  in  the  equa¬ 
tions  for  yvz  and  yv±,  and  the  drift  due  to  the  gradient  of  the  external  field 
in  the  equations  for  x  and  y.  The  latter  equations  do  not  include  the  curva¬ 
ture  drift  but  contain  additional  terms,  due  to  the  fact  that  the  axis  z  is 
defined  to  be  the  axis  of  the  cavity  and  not  the  exact  direction  of  the  exter¬ 
nal  field. 


Similarly,  the  slow-time-scale  equations  for  the  growth  rate  and 
frequency  shift  of  each  eigenmode  are: 
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(12) 


x<^Ji  <+  l)(^  cosC±  -  knSnCm  Sinc±)^>  , 

where  <5nm  Q  =  1,  if  either  n  =  0  or  m  =  0,  and  6nm  Q  =  0,  otherwise,  and  the 

angular  brackets  denote  an  average  over  all  particles  inside  the  cavity,  their 

number,  i'  v,  being  accounted  for  by  the  beam  plasma  frequency,  co„. 

max  p 

Let  us  review  briefly  the  steps  involved  in  the  derivation  of  Eq.  (12). 


The  wave  equation  was  applied  on  the  cavity  fields,  expressed  by  the  summa¬ 
tions  in  Eq.  (3).  The  orthogonality  of  the  eigenmodes  was  used  to  decouple  the 
wave  equation,  by  integrating  over  the  cavity  volume  after  multiplying  by  the 


appropriate  tri gonometri c  functions,  with  the  term  1  +  6nm  n  arising  out  of 
the  different  normalization  when  either  n  =  0  or  m  =  0.  The  boundary  condition 
corresponding  to  the  power  radiated  outside  the  cavity  is  accounted  for  by 
the  quality  factor,  Qnm£»  which  depends  on  the  cavity  and  the  structure  (but 
not  the  amplitude)  of  each  eigenmode.  The  terms  r/Q  +  (r2  +  dr/dt)/w  and 
dA/dt  are  neglected  compared  to  A  and  to  w/Q  +  2r,  respectively,  and  a  multi¬ 
plication  by  cos<)>  and  sin<}>  with  omission  of  the  fast  terms  leads  to  Eq.  (12). 
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5.  ANALYTICAL  APPLICATION  IN  A  WAVEGUIDE 

Before  using  Eqs.  (10)  and  (12)  for  any  extensive  study  either  to  predict 
the  behavior  of  a  given  oscillator  or  to  investigate  the  physical  processes 
involved  in  multi -mode  interactions,  it  is  prudent  to  first  establish  relia¬ 
bility  in  their  correctness.  This  purpose  is  served  by  the  present  section, 
in  which  our  equations  are  used  to  describe  analytically  a  variety  of  situations, 
in  order  to  compare  these  results  with  the  appropriate  earlier  ones.  The  only 
configuration  in  our  rectangular  geometry,  for  which  sufficient  information 
exists  to  compare  to  our  results,  is  that  of  an  amplifier  operating  at  a 
single  mode.  The  conversion  of  our  equations  to  that  configuration  is  straight¬ 
forward:  it  suffices  to  eliminate  t_,  the  backward  wave  dependence,  to  remove 
the  meaningless  term  oj/2Q  in  the  equation  for  r  and  to  divide  the  equations 
for  r  and  A  by  2,  since  their  derivation  for  the  amplifier  case  involves  an 
integration  only  over  the  cross  section  of  the  waveguide,  and  k£  has  now  an 
arbitrary  value,  not  quantized  by  the  meaningless  zmax. 

First  it  is  noted  that  the  equations  for  the  electrons,  Eq.  (10),  for  the 
special  case  of  no  nonuniformity  reproduce  earlier  results  (3),  when  m=0, 
n  f  0,  and,  with  an  appropriate  rotation  of  the  coordinate  axes,  when  n=0, 
m  f  0.  On  the  other  hand,  as  noted  in  the  previous  section,  the  effects 
associated  with  the  nonuniform  external  field,  that  is,  the  acceleration  due 
to  a  longitudinal  gradient,  the  drift  due  to  a  transverse  gradient  and  the 
drift  due  to  an  imperfect  alignment  of  the  external  field  with  the  cavity 
(or  waveguide)  axis,  are  properly  accounted  for  by  our  equations. 

In  order  to  see  whether  Eq.  (.12)  describes  the  correct  evolution  of  the 
fields,  and  because  of  the  usefulness  of  a  dispersion  relation  for  arbitrary 
m  and  n.  It  is  advisable  to  derive  such  a  general  relation  and  to  compare  it 
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to  the  earlier  result  (2)  for  m=0.  In  order  to  obtain  the  dispersion  relation, 
it  is  assumed  that  r  and  A  are  constant.  The  equations  of  motion  are  solved 
self-consistently  with  the  exponentially  amplified  field  and,  as  appropriate 
for  the  high  gain  regime  (It  »  1),  the  small  terms  are  dropped.  A  final 
integration  over  the  initial  phase  angle  of  the  electrons,  assumed  uniformly 
distributed,  gives  the  expressions  for  r  and  A,  which  are  combined  to  yield 
the  dispersion  relation, 


a2  -  (k  „c)2  - 
'  nmfc  ' 
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where  fQ(x,  y,  pz,  px)  is  the  unperturbed  distribution  function  and 
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K2  =  (k  C  S  )2  +  (k  S  C  )2  . 

'mnnr  1  n  n  nr 

Introducing  into  Eq.  (13)  a  monoenergetic  pencil  beam,  f  =  (*max  ymax)/(2lTPx) 

«(x  "  *max/2)  "  ^max/2)  6(Px  “  PiJ  6(p|j  "  p||  oJ  reProduces  earlier 

expressions.  Furthermore,  omitting  the  less  important  second  term  in  brackets 
in  Eq.  (13)  and  confining  our  attention  to  the  behavior  of  the  dispersion  rela¬ 
tion  away  from  the  isolated  (and  uninteresting)  real  root  a  =»  -  k^c,  gives 

the  simplified  cubic  relation 
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where  0+q  =  ft2/y0  -  knmJl  +  k£vz  is  the  frequency  mismatch  and  the  coupling 
parameter,  aQ,  is  given  by 
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As  indicated  by  Eq.  (15),  for  parameters  sufficiently  above  threshold  and  for 

adequate  synchronism,  the  evolution  of  the  instability  depends  on  the  parameters 

0,  and  of. 

+o  o 

An  important  approximation  in  the  derivation  of  Eq.  (13),  the  dispersion 
relation,  has  been  that  the  small  argument  expansion,  J-j(A)  =  A/2,  has  been 
employed,  where 
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This  approximation  is  appropriate  for  the  linear  regime  before  saturation  is 
reached,  i.e.,  when  flnmj(and  u£)  is  small.  Since  this  approximation  becomes 
increasingly  inaccurate  after  A  =  1.84  (where  d^/dA  =  0)  and  definitely  fails 
when  A  =  3.83  (where  J-j  =0),  these  two  values  give  the  range  of  A  in  which 
saturation  is  expected,  1.84  <  A  <  3.83,  and  therefore  also  an  estimate  for 
the  saturation  amplitude,  ft  ^  from  Eqs.  (17)  and  (18). 
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A  final  application  of  our  equations  is  for  the  case  of  a  saturated 
monochromatic  wave,  with  almost  constant  amplitude  and  frequency.  Assuming 
that  va  and  y  vary  only  weakly  (low  efficiency),  then  the  most  important 
variation  is  that  of  the  phase  angle  c+,  which  is  found  to  satisfy  the  pendulum 
equation, 

d2? 

jp-  =  sin(s+  +  n)  .  (19) 


where  the  bounce  frequency,  has  been  defined  in  Eq.  (18),  and  the  reference 
angle  n  is  defined  by  the  field  structure. 
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Both  oob  and  n  simplify  to  the  known  expressions  (3)  when  m=0,  n^O.  Furthermore, 
the  corresponding  equations  for  the  growth  rate  and  frequency  shift  of  the 
wave  give 
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In  view  of  the  electron  bunching  that  has  occurred  during  the  exponential 
stage  of  the  evolution,  the  evolution  of  the  center  of  mass  of  the  bunch 
during  the  nonlinear  stage  according  to  the  pendulum  equation,  Eq.  (19), 
schematically  represented  in  Fig.  3,  gives  from  Eq.  (21)  the  characteristic 
one-  and  two-hop  variation  of  r  and  A  during  one  bounce  period,  shown  in  Fig.  4. 
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NUMERICAL  APPLICATION  IN  A  WAVEGUIDE 

In  the  previous  section  we  have  derived  equations  governing  various 
aspects  of  monochromatic  wave  evolution  in  a  waveguide.  We  have  shown  there 
that  these  equations  reproduce  earlier  results  in  the  appropriate  limit  ufO, 
n^O,  thus  establishing  by  implication  the  correctness  of  Eqs.  (10)  and  (12). 

The  next  subject,  covered  by  this  section,  is  to  establish  that  the  algorithm 
based  on  these  equations  gives  acceptably  accurate  numerical  results.  In  all 
simulations  discussed  in  this  section  a  low  amplitude  monochromatic  wave  was 
initialized  in  a  system  with  periodicity  length  zmgY  *  The  initial 

distribution  of  the  electron  position  in  z  and  the  phase  angle  i|>  was  uniform, 
with  the  initial  values  assigned  to  each  electron  by  the  use  of  a  random  number 
generator,  to  add  an  experimental  flavor  to  the  results,  while  x,  y,  v  and  vz 
were  initially  the  same  for  all  particles.  The  subsequent  evolution  of  the 
electron  and  wave  properties  was  obtained  by  a  successive  application  of  the 
long  time-scale  algorithm  of  Eqs.  (10)  and  (12). 

The  time  evolution  of  the  amplitude,  growth  rate  and  frequency  shift  is 

shown  in  Fig.  5,  for  a  simulation  with  kp  =  0.958  Oz/c  (n  =  1,  xmax  =  3.28  c/«z), 

km  =  °»  k£  =  °-108  Vc  tzmax  =  58,0  c/nz^  “p  =  10~3  °z»  vzo  =  0,22  c’ 

v±Q  =  0.26  c,  y0  =  1.0636  and  xQ  =  0.5  xpiax.  For  these  values,  the  theoretical 
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predictions  are  T  =  1.68  x  10  nz,  A  a  0.97  x  10  Q  ,  to  which  the  wave 

evolution  is  in  good  agreement  in  the  linear  regime.  The  estimate  Amax  *  1.84 

-4  -4 

gives  «max  =  1.16  x  10  ,  comparable  to  the  measured  value  =  1.30  x  10  . 

Finally,  at  saturation,  oscillations  are  observed  at  the  bounce  frequency 

Wjj  -  2.72  x  10  «z,  as  expected  from  Eq.  (.21).  Some  skewness  in  the  variation 

of  r  and  the  inequality  of  the  peaks  of  A  in  Fig.  5b,  are  accounted  for  by  the 

dependence  of  Eq.  (21)  on  the  inverse  of  the  wave  amplitude,  which  was  assumed 

to  be  constant  when  Fig.  4  was  prepared. 
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ig.  5  -  Evolution  of  (a)  the  amplitude  and  (b)  the  growth  rate  and  the 
frequency  shift  for  the  reference  simulation. 


Some  additional  simulations  are  presented  in  Table  1.  Case  0  is  the 
simulation  already  described,  which  we  regard  as  the  reference  simulation. 

The  simulations  of  Table  1  have  the  same  initial  values  as  the  reference 
case,  except  as  indicated  in  the  input  column  and  discussed  below.  As  it 
is  seen,  cases  1-5  have  the  appropriate  scaling  with  regard  to  up  and  v^. 

The  absence  of  amplification  in  case  5  is  consistent  with  the  value  of  v 
being  below  threshold.  In  cases  6  and  7,  where  the  distribution  of  the 
guiding  center  is  taken  as  f(x)  =  1/tt  and  (sinx)/2  respectively,  rather  than 
xmax  ”  xmax'^’  the  9rowtf)  rates  and  frequency  shifts  are  reduced,  since 
the  average  electron  is  acted  upon  by  a  weaker  field.  No  amplification  occurs 
when  a  delta  distribution  is  introduced  at  the  node  of  the  n=2  field  (case  8), 
while  case  9  gives  essentially  the  same  results  as  case  6,  consistent  with  the 
fact  that  the  average  field  is  independent  of  the  number  of  nodes  for  a  uni¬ 
form  particle  distribution  in  x. 

The  last  set  of  numerical  results  are  summarized  in  Fig.  6.  In  these 
simulations  the  purpose  was  to  check  the  extent  to  which  the  dependence  of  r 
and  A  and  (at  saturation)  satisfy  the  dispersion  relation  of  Eq.  (15) 
and  the  estimate  A  =  1.83,  for  various  values  of  the  frequency  mismatch,  e+. 
The  data  in  Fig.  6  were  obtained  with  input  parameters  equal  to  those  of  the 
reference  simulation,  except  for  the  values  of  Jlz  or  the  periodicity  length, 
in  order  to  produce  a  variety  of  values  for  e+.  The  quantities  in  Fig.  6  are 
scaled  to  the  relevant  parameter,  aQ,  defined  in  Eq.  (16).  As  can  be  seen, 
given  the  difficulty  of  measuring  these  values  from  curves  like  those  of 
Fig.  6,  the  agreement  is  very  good  with  regard  to  r  and  A,  and  the  saturation 
amplitude  only  slightly  exceeds  the  value  given  by  the  estimate  A  =  1.83. 
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Table  1 


DEPENDENCE  OF  r,  A  and  nmax  ON  THE  INPUT  PARAMETERS 
FOR  APPLICATIONS  IN  A  WAVEGUIDE 


Case 

Input 

r/ftz 

A/nz 

n  /n 

max'  z 

0 

Reference 

1.7  x 

10'3 

1.0  x 

10‘3 

1.3  x  10“4 

1 

x,e 

so 

N 

II 

O 

1 

4.0  x 

10'4 

2.4  x 

10‘4 

6.7  x  10"6 

2 

cytoz  =  10‘2 

7.7  x 

10-3 

4.5  x 

IQ'3 

3.2  x  10‘3 

3 

yjc  =  0.6 

2.1  x 

10"3 

1.4  x 

10"3 

2.3  x  10"4 

4 

< 

h 

O 

It 

o 

8.2  x 

10‘4 

5.0  x 

10’4 

1.0  x  10‘4 

5 

vjc  =  0.01 

0 

0 

undef. 

6 

uniform  distr. 

1.2  x 

10'3 

5.7  x 

10-4 

1.3  x  10‘4 

7 

sine  distr. 

1.6  x 

10-3 

6.2  x 

10"4 

1.4  x  10'4 

8 

n  =  2,  delta  distr. 

0 

0 

undef. 

9 

n  =  2,  uniform  distr. 

1.2  x 

10-3 

5.3  x 

10"4 

1.2  x  10'4 
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7.  EFFECTS  OF  AN  EXTERNAL  FIELD  NONUNIFORMITY 

It  has  been  shown  theoretically  that  the  efficiency  enhancement  is  achieved 
if  the  external  field  has  some  nonuniformities,  either  in  the  form  of  a 
step(15)  or  of  a  constant  gradient  (16).  The  same  effect  has  been  also  shown 
experimentally  (17).  A  qualitative  explanation  of  this  phenomenon  is  that  the 
introduction  of  the  nonuniformity  keeps  the  frequency  mismatch,  0+,  at 
favorable  values  over  a  prolonged  interval,  thus  delaying  the  occurrence  of 
saturation.  Some  theoretical  considerations  of  the  effects  associated  with 
the  nonuniformity  are  given  in  this  section.  Although  these  considerations 
are  of  introductory  nature,  they  may  prove  useful  for  the  elucidation  of  the 
interaction. 

In  order  to  simplify  the  analysis,  we  consider  the  simple  case  in  which 
the  transverse  velocity  is  not  too  much  perturbed  during  the  interaction, 
e.g.,  during  the  transit  time  in  the  cavity.  This  assumption,  already  employed 
in  Sec.  5  for  the  derivation  of  Eq.  (.19),  is  rigorously  applicable  only  when 
the  efficiency  of  conversion  is  low,  however,  it  is  felt  that  the  results  may 
be  extrapolatable  to  the  high  efficiency  regime,  at  least  qualitatively.  In 
the  sense  of  this  approximation,  the  equation  for  c+  is  now  given  by 


d25+ 

-  2 
dt*=  wb 


(22) 


where  is  defined  in  Eq.  (18)  and  the  effects  of  the  nonuniformity  are 
described  by  the  nonuniformity  ratio,  R,  expressing  the  gradient  of  the 
external  field  in  relation  to  the  wave  field  and  defined  by 


I 


Generally,  the  difference  of  the  last  two  terms  in  the  above  definition  is 
small  compared  to  the  first  term,  so  that  approximately  Ro>£  =»  <S(Qz/y)/6t, 
where  6t  is  the  interaction  duration  and  <$(Qz/y)  is  the  corresponding  change 
in  the  electron  gyrof requency . 

Equation  22  represents  the  pendulum  equation,  with  an  added  external 
torque  expressed  by  R.  Several  features  associated  with  this  equation  are 
presented  in  Ref.  22.  Here  it  suffices  to  note  that  the  condition  for  existence 
of  trapped  particles  is 

|R|  <  1  .  (24) 

and  that  the  separatrix  extends  approximately  to  the  values 

1 9+1  -  2  a»b(l  -  |Rj)  .  ‘  (25) 

For  the  representative  case  R  =-0.4,  the  trajectories  of  the  particles  are 
shown  in  Fig.  7b,  in  phase  space  coordinates  c+  +  n  and  e+  =  dc+/dt.  These 
trajectories  are  to  be  compared  to  the  case  of  a  uniform  external  field. 

Fig.  7a,  where  by  definition  R  =  0.  In  both  cases,  trapped  particles  like 
1  and  V  execute  oscillatory  orbits  about  the  vortex,  which  is  located  at 
C+  +  n  =  tr,  when  R  =  0,  and  at  t+  +  n  =  ir  +  arcsin  R,  when  R  i  0.  On  the 
other  hand,  untrapped  particles  like  2  and  2'  execute  a  wavy  motion,  when 
R  =  0,  never  crossing  the  resonance  line  0+  *  0,  while  the  introduction  of  a 
nonzero  value  for  R  introduces  an  additional  upward  drift,  causing  untrapped 
particles  to  eventually  cross  the  resonance  line. 

Based  on  the  above  remarks,  but  before  describing  the  mechanism  responsible 
for  the  improved  efficiency  when  R  f  0,  it  is  convenient  to  discuss  first  the 
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uniform  case,  R  =  0,  within  the  framework  of  our  equations.  Recalling  from 
Eq.  (21)  that  energy  transfer  from  the  particles  to  the  wave  is  associated 
with  <  sin(s+  +  n))  >  0,  that  is,  when  the  average  phase  is  in  the  interval 
0  <  C+  +  n  <  tt,  we  note  from  the  direction  of  the  trapped  particle  motion 
in  Fig.  7a  that  this  interval  is  satisfied  for  interaction  times,  t,  such 
that  wbT  <_tt,  provided  that  0^  <  0.  This  is  also  consistent  with  the  fact 
that  during  such  a  time  interval,  the  value  of  0+  =  Qz/y  -  w  +  k^vz  increases 
for  these  particles,  resulting  in  a  decrease  of  y  and  energy  loss.  The 
untrapped  however,  do  not  contribute  to  the  average  of  sin(c+  +  n)  substantially, 
as  can  also  be  seen  from  the  fact  that  their  value  of  0t  suffers  weaker  changes 
compared  to  the  trapped  particles.  Accordingly,  the  optimal  situation  is 
determined  by  a  trade-off  between  the  desire  for  a  large  value  for  -  0^,  so 
that  the  average  of  sin(c+  +  n)  is  large  for  the  trapped  particles,  and  for 
a  small  value  for  -  0+Q,  in  order  to  reduce  the  percentage  of  idle  untrapped 
particles. 

Let  us  now  turn  our  attention  to  the  additional  effects  caused  by  the 
nonuniformity.  We  consider  first  the  case  R  >  0  (but  R  <  1),  that  is  of  an 
external  field  increasing  in  the  direction  of  the  motion.  This  case  is 
characterized  by  the  transition  of  the  untrapped  particles  through  the  line 
0+  =  0  preferentially  at  values  satisfying  sin(c+  +  n)  >  0,  while  the 
trapped  particles  oscillate  about  the  vortex,  where  sin(c+  +  n)  <0. 

Accordingly,  the  desirable  initial  properties  of  the  beam  are  those  that 
correspond  to  untrapped  particles.  On  the  other  hand,  distant  untrapped 
particles  are  not  particularly  attractive,  since  such  particles  will  first 
execute  a  number  of  oscillations  in  c+  +  n  =  (0,2tt),  before  crossing  the 
resonance  line,  0+  =  0,  with  sin(c+  +  n)  >  0.  Therefore,  the  most  effective 
initial  arrangement  for  the  beam  is  that  for  which  0^=  -  2^  (1-R),  i.e.. 


the  beam  is  initially  close  to  the  edge  of  the  trapping  region.  Following 
the  trajectories  of  the  particles  for  such  an  initial  arrangement,  e.g., 
according  to  Eq.  (22)  or  Fig.  7b,  one  can  see  that  a  bunch  is  formed,  which 
remains  over  extended  times  in  the  region  sin(c+  +  n)  >  0.  The  effect  of  this 
bunch  can  be  optimized  by  adjusting  the  value  of  R  and  e+Q.  However,  it  can 
be  seen  that  the  corresponding  optimized  efficiency  will  be  larger  when  R  +  0 
than  when  R  *  0,  since  the  bunch  is  formed  out  of  all  particles,  and  not  only 
out  of  that  fraction,  which  is  initially  trapped,  as  is  the  case  when  R  =  0. 

The  opposite  situation,  that  of  an  external  field  gradient  opposite  to 
the  direction  of  motion  (R  <  0,  but  R  >  -  1),  is  characterized  by  a  location 
of  the  vortex  such  that  sin(c+  +  n)  >  0.  Therefore,  in  this  case  it  is 
advantageous  to  arrange  for  a  bunch  to  be  formed  out  of  the  trapped  particles, 
as  shown  in  Fig.  8b.  Comparing  this  evolution  to  the  uniform  case,  it  can 
be  seen  that  this  bunch  is  expected  to  be  weaker  than  that  of  the  uniform  case, 
since  the  trapping  region  is  smaller  and  relatively  more  particles  will  be 
untrapped  and  therefore  idle,  as  far  as  the  energy  conversion  is  concerned. 

On  the  other  hand,  this  bunch,  although  weaker,  will  be  located  closer  to  the 
convenient  phase  angle,  c+  +  n  »  ir  +  arc  sin  R,  when  R  <  0.  If  this  shift  of 
the  vortex  (and  bunch)  phase  is  stronger  than  the  reduction  of  the  bunch 
density,  then  the  efficiency  for  R  <  0  will  be  larger  than  the  uniform  case. 
However,  this  question  cannot  be  answered  within  the  present  qualitative 
discussion.  It  is  clear,  nevertheless,  that  if  the  trapped  particles  are 
allowed  to  execute  a  number  of  bounce  oscillations,  then  they  will  contribute 
energy  to  the  wave  during  each  oscillation,  since  for  the  vortex  sin(c+  +  n)  * 

-  R  >  0,  hence  an  improved  efficiency  is  expected  when  »  n,  i.e.,  for 
systems  with  high  field  amplitudes  (e.g.,  driven  by  high  currents),  or  long 
transit  times,  t,  e.g.,  an  amplifier. 
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the  bunch  for 


j  It  is  interesting  to  note  that  when  R  f  0,  limitations  associated  with 

beam  quality  are  not  so  severe  as  when  R  =  0.  In  the  uniform  case,  a  beam 

j  temperature  associated  with  a  spread  of  6+ ,  such  that  66+  >  essentially 

* 

kills  the  instability,  since  in  this  case  the  effect  of  a  particle  with 

j  initial  values  0+o,  c+£)  +  n  is  compensated  by  that  of  a  particle  at  -  0+o, 

r 

|  2tt  -  (s+Q  +  n)-  This  is  not  the  case,  however,  when  a  nonuniformity  is  intro- 

[  duced.  For  example,  when  R  >  0,  the  external  field  can  be  adjusted  so  that 

t 

|  at  the  entrance  of  the  cavity  0+  <  -  2oib  (1  -  R)  for  all  particles.  The 

j  ensuing  motion  of  these  untrapped  particles  (Fig.  9a)  will  bring  them  near 

f 

resonance,  where  they  will  occupy  a  region  with  <sin(c+  +  n)>  >0.  transferring 
energy  to  the  wave.  Similarly,  when  R  <  0,  the  initial  value  of  0+  can  be 

I  arranged  to  be  on  the  average  close  to  zero.  Then,  the  untrapped  particles 

i  1 

,  |  will  eventually  leave  the  resonance  region,  spread  uniformly  in  c+,  while  the 

|  trapped  particles  will  remain  inside  the  separatrix,  with  <sin(c+  +  n)>  >  0. 


|  The  case  R  >  0  is  particularly  adoptable  to  a  cavity,  since  enrgy  conversion 

j 

1  is  limited  to  the  time  it  takes  the  untrapped  particles  to  move  around  the 

separatrix,  while  the  case  R  <  0  is  primarily  applicable  to  an  amplifier,  to 

i 

take  advantage  of  the  extended  time,  over  which  the  trapped  particles  are  able 
to  convey  energy  to  the  wave.  In  either  case,  calculations  for  the  correspond¬ 
ing  efficiencies  can  be  made  based  on  the  data  of  Ref.  22. 

( 

Finally,  it  can  be  seen  that  mathematically  the  introduction  of  a  gradient 
in  the  external  field  is  equivalent  to  a  variation  of  the  cavity  cross  section, 
1  since  both  cause  0+  to  vary  independently  of  the  cavity  fields.  A  generalized 


4  expression  can  be  written  for  R, 


Schematic  representation  for  the  evolution  of  a  beam  with  temperature 
(a)  when  R  >  0,  with  no  particles  initially  trapped,  and  (b)  when 
R  <  0,  with  the  trapping  region  initially  filled. 


with  the  derivative  being  evaluated  in  the  absence  of  the  cavity  fields.  When 
dnz/dz?  0,  Eq.  (26)  yields  the  expression  already  presented  (Eq.  (23)).  A 
variable  cross  section  causes  to  be  a  function  of  z,  since  k|  a  (w/c)2  - 
"2(n2/xLw(z)  +  m2/yLv(z))»  hence  the  corresponding  expression  for  R  is 

iTlaX  maX 


(27) 


In  addition,  an  alternate  situation  where  Eq.  (27)  is  applicable  is  when  the 
cavity  is  lined  with  a  thin  and  weak  dielectric  material,  with  its  dielectric 
constant  varying  with  z,  while  an  appropriate  expression  can  be  obtained  from 
Eq.  (26)  when  a  weak  electrostatic  field  is  externally  applied  along  the  axis 
of  the  cavity. 
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8.  APPLICATION  TO  AN  OVERMODED  CAVITY 


In  this  section  we  present  and  discuss  various  results  obtained  from  the 

numerical  application  of  the  algorithm  formed  from  Eqs.  (10)  and  (12)  to  the 

case  of  a  rectangular  cavity,  in  which  more  than  one  eigenmode  is  in  close 

resonance  with  the  beam.  We  have  chosen  the  parameters  so  that  they  correspond 

as  closely  as  possible  to  those  of  typical  NRL  experimental  arrangements.  Thus, 

the  electron  beam  was  chosen  with  initial  velocities  vZQ  =  0.267c  and  v^  ■ 

0.400  c,  giving  a  relativistic  factor  of  y0  55  1.14  and  an  anisotropy  ratio 

vjx/vzo  =  typical  of  the  NRL  beam.  The  dimensions  of  the  rectangular 

cavity  were  zmav  =  33.10  c/O,  in  length,  xma„  =  6.78  c/O,  in  width,  with  the 

height,  yma  .  sufficiently  short,  so  that  only  the  modes  with  m  =  0  can  be 

resonant.  In  this  cavity,  resonant  modes  are  those  with  n  =  2,  hence 

kn  *  0.927  Oz/c,  and  with  £  =  1  and  £  =  2,  corresponding  to  k&  =  £  x  0.095  Oz/c. 

This  choice  gives  a  ratio  of  kn/k£  equal  to  the  corresponding  ratio  in  the  NRL 

cylindrical  cavity  with  a  length  to  diameter  ratio  of  4,  with  the  coordinates 

(x,  y,  z)  corresponding  to  (r,  6,  z)  of  the  cylindrical  geometry  and  xmflX 

corresponding  to  the  diameter.  The  beam  is  placed  at  x  =  0.33  xm_„,  since  at 

that  point  the  ratio  E^/Bz  =  tan(knx)  =  1.84  is  equal  to  the  corresponding 

value  at  EQ/BZ  for  the  cylindrical  cavity  at  the  point  where  Efi  is  maximum. 

The  corresponding  values  of  the  empty  cavity  frequencies  are  =  0.932  «z  and 

=  0.946  n2  ,  the  frequency  mismatches  relative  to  the  foreward  components  of  these 

-2  -2 

modes  are  6+1  *  -  2.95  x  10  flz  and  0+2  =  -  1.85  x  10  «z»  where  the  subscripts 
refer  to  £,  the  axial  harmonic  number.  The  transit  time  of  the  electron 
through  the  cavity  is  t  ■  zmax/v20  =  124  and  at  the  position  of  the  beam 
the  effective  transverse  wavenumber,  defined  by  Eq.  (14c),  is  K  *  kR  s1n(knx) 

*  0.811  nz/c. 

For  the  numerical  Implementation  of  our  slow- time-scale  algorithm,  the 
cavity  is  Initially  assumed  empty,  with  low  amplitude  standing  cavity  waves. 
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The  beam  electrons  are  Initialized  all  with  the  same  values  for  x,  vz,  yA,  as 
given  above,  while  their  initial  position,  zQ  <  0,  and  phase  angle,  t|>0,  are 
given  a  uniform  distribution  by  use  of  the  random  number  generator.  The 
electrons  advance  according  to  Eq.  (10),  if  Inside  the  region  of  interaction, 

0  <  z  <  zmaw,  and  according  to  unperturbed  cyclotron  motion,  given  by  Eq.  (6) 
when  outside  the  cavity.  Depending  on  the  density  (or  current)  that  the 
simulation  electrons  represent,  it  was  found  that  approximately  100  to  300 
electrons  suffice  to  yield  accurate  statistics,  while  an  accurate  integration 
of  their  equation  of  motion  is  obtained  with  a  time  step.  At,  ranging  from 
1/8  to  1/15  of  the  transit  time.  Because  of  the  continuous  positioning  of 
the  electrons  along  the  z-axis,  they  do  not,  in  general,  cross  the  boundaries 
at  the  quantized  multiples  of  At,  hence  for  electrons  near  the  boundaries  the 
algorithm  is  applied  over  the  effective  time  of  interaction,  that  is,  the 
fraction  of  At  which  corresponds  to  motion  inside  the  cavity.  Finally,  since 
electrons  that  cross  the  exit  boundary  do  not  contribute  to  the  interaction 
any  more,  they  are  "destroyed"  and  their  memory  locations  in  the  computer  is 
occupied  by  properties  of  electrons,  which  are  about  to  enter  the  cavity.  This 
procedure  results  in  substantial  savings  on  the  size  of  memory  required,  since 
the  transit  time  is  a  very  small  fraction  of  the  total  time  required  to  observe 
a  semblance  of  a  steady  state.  In  implementing  this  procedure,  care  has  been 
exercised  to  assure  that  a  uniform  particle  density  in  z  is  maintained,  before 
they  enter  the  cavity. 

In  Table  2,  we  present  the  results  of  ten  simulations  for  the  cavity  and 
modes  presented  above.  In  these  simulations  the  quality  factors  are  Q-j  *  1000 
and  Q2  =  700,  and  the  simulation  electrons  represent  different  plasma  densities, 
as  given  by  the  value  of  up.  In  cases  1-6  and  10,  both  modes  are  Initialized 
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Table  2 


\  0 


4.71  5.58  1.18  6.91  0.53  14.6 


with  a  low  amplitude,  ft  =  10  ftz>  and  one  mode  became  dominant,  as  denoted 
In  the  Table.  This  mode  was  initially  amplified  essentially  exponentially, 
at  a  rate  given  by  r,  and  later  approached  an  oscillatory  asymptotic  satura¬ 
tion,  with  the  value  ftmai/  and  the  corresponding  value  for  the  bounce  frequency, 
as  given  in  the  Table.  Additional  entries  in  the  Table  are  the  ratio 
ft^/oip,  the  product  of  the  bounce  frequency  times  the  transit  time,  the 
ratio  -  6+/t0jj  of  the  frequency  mismatch  to  the  bounce  frequency,  and  the  effi¬ 
ciency,  n,  of  conversion,  defined  by  Eq.  (4)  as  the  ratio  of  radiation  to  beam 
entry  power.  Caser  7-9  are  a  repetition  of  cases  10,  1,  and  5,  except  that 

the  nondominant  mode  of  the  original  cases  was  initially  primed  to  a  value  of 
_2 

10  ftz  for  the  amplitude.  This  resulted  in  these  modes  becoming  artificially 
dominant,  with  parameters  at  saturation  as  given  in  the  Table. 

Typically,  the  evolution  of  the  amplitude  and  frequency  of  the  two  modes 
is  as  shown  in  Fig.  10,  where  the  actual  values  of  the  parameters  are  those  of 
case  4.  The  dominant  mode,  which  here  is  that  with  Z= 2,  is  seen  to  grow 
exponentially,  with  a  weak  oscillation  superposed,  and  later  to  approach  an 
asymptotic  state.  The  nondominant  mode  initially  attempts  to  grow,  at  a 
weaker  rate  and  with  relatively  stronger  oscillations,  reaching  a  maximum 
amplitude  with  a  value  depending  on  the  beam  density.  The  peak  amplitude  is 
always  reached  near  the  time  of  transition  from  the  exponential  amplification 
to  approach  to  equilibrium  of  the  dominant  mode.  Beyond  that  time,  the  non¬ 
dominant  mode  decays  on  the  average,  while  the  oscillations  persist.  In  all 
simulations  the  period  of  these  oscillation  is  equal  to  450  ft"1,  equal  to 
the  beat  period,  2ir/(a>2  -  oj-j ) ,  of  the  two  modes.  Oscillations  of  the  same 
period  are  also  present  in  the  amplitude  of  the  dominant  mode,  however  these 
oscillations  are  too  weak  to  be  seen  in  Fig.  10.  With  regard  to  the  frequency 
plots,  one  first  notes  a  relatively  wide  dispersion  in  the  data  points. 


Fig.  10  -  Evolution  of  the  (a)  amplitude  and  (b)  frequency  of  the  resonant 

modes  for  case  4  of  Table  2.  The  horizontal  broken  lines  In  (b) 
give  the  empty  cavity  frequencies. 


primarily  during  the  initial  evolution  of  the  i,«*2  mode  and  always  for  the 
£=1  mode.  The  dispersion  is  directly  related  to  the  fluctuations  associated 
with  the  randomness  in  the  particle  initialization.  Similar  fluctuations  are 
present  in  the  evolution  of  the  growth  rates.  Such  fluctuations  are  of  no 
major  concern,  since  they  affect  only  the  instantaneous  values,  r  and  A, 
while  particle  dynamics  are  governed  primarily  by  the  wave  amplitudes  and 
phases,  which  are  much  smoother  (see  Figs.  10a  and  11)  being  integrated  values. 
Apart  from  these  fluctuations,  it  can  be  seen  that  the  frequency  of  the 
dominant  mode  essentially  approaches  a  constant  value,  shifted  from  the  empty 
cavity  value  by  approximately  1%,  while  the  frequency  of  the  £=1  mode  has 
similar  beat  oscillations,  like  the  amplitude. 

In  all  cases  of  Table  2,  and  in  many  more  two-mode  cases  not  presented 
here,  the  asymptotic  behavior  of  the  spectrum  is  dominated  by  one  mode,  the 
amplitude  and  energy  of  the  additional  mode  being  respectively  one  and  two 
orders  of  magnitude  smaller.  In  such  cases,  it  is  appropriate  to  treat  the 
nondominant  mode  as  a  stable  perturbation  driven  by  a  constant  amplitude 
monochromatic  field,  while  the  dynamics  of  the  electrons  is  determined  by  the 
large  amplitude  field  only.  If  we  assume  that  the  motion  of  the  particles  is 
adequately  represented  by  the  pendulum  equation,  Eq.  (19),  in  the  limit  of 
small  oscillations  about  +  n  =  *,  then  integrating  sin(t+  +  n)  over  the 
transit  time,  substituting  this  time  average  in  place  of  the  ensemble  average 
In  Eq.  (21a),  with  r  replaced  by  w/2Q,  as  is  appropriate  for  a  steady  state 
oscillator,  and  using  also  the  definition  of  u^,  Eq.  (18),  with  =  0.5, 
one  obtains  a  relation  for  the  amplitude  of  the  dominant  mode. 


and  accordingly,  the  efficiency  Is  obtained  from  Eq.  (4), 
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The  above  equations  demonstrate  the  known  fact  that  the  wave  frequency  must 
exceed  the  gyrofrequency,  i.e.,  6+0  <  0,  in  order  for  positive  energy  conver¬ 
sion  to  be  possible.  The  coefficient  f  has  the  value  f  =  1  in  the  idealized 
situation  when  the  pendulum  equation  is  applicable,  the  bunch  is  perfect  and 
consists  of  all  particles,  and  the  bunch  performs  half  of  a  bounce  oscilla¬ 
tion  during  transit.  Deviations  from  such  idealized  conditions  reduce  somewhat 
the  value  of  f.  For  the  simulations  already  discussed,  the  idealized  values 
of  from  Eq*  (28)  when  f  =  1,  are  2.24  and  1.45,  for  the  £=1  modes 

respectively,  and  as  can  be  seen  from  Table  2,  depending  on  and  6+/^,  the 
actual  values  of  ft^/Wp  indicate  that  f  £0.6.  The  same  holds  true  also  for 
the  simulations  not  presented  here. 


Let  us  turn  our  attention  now  on  the  effects  which  the  large  amplitude 
wave  has  on  the  nondominant  one.  Since  the  bunch  is  attached  to  the  dominant 
mode,  the  current  it  represents  will  altematingly  act  constructively  or  des¬ 
tructively  on  the  secondary  mode,  depending  on  whether  they  are  in  or  out  of 
phase.  This  succession  is  repeated  at  the  beat  period,  as  observed  in  the 
simulations.  Similarly,  due  to  the  oscillations  of  the  reactive  component  of 
the  current  relative  to  the  low  amplitude  wave,  its  frequency  oscillates  about 
the  empty  cavity  value.  This  can  be  seen  in  Fig.  10b  up  to  the  time  t  *  3000. 
However,  beyond  that  time,  which  incidentally  coincides  with  the  time  the 
dominant  mode  amplitude  has  practically  reached  the  final  saturation  value, 
the  frequency  of  the  driven  mode  appears  to  oscillate  about  that  of  the 
dominant  mode.  This  can  be  better  seen  in  Fig.  11,  where  the  phase  of  mode 
t=l  relative  to  that  of  S.=2,  -  <J>2»  is  plotted  against  time,  the  latter  in 
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an  expanded  scale  compared  to  Fig.  10,  for  improved  clarity.  It  can  be  seen 
there  that  the  relative  phase  initially  executes  a  number  of  complete  rota¬ 
tions,  during  which  it  preferentially  remains  a  relatively  longer  time  in  the 
vicinity  of  -  <j>2  =*  it,  at  times  which  coincide  with  the  occurrence  of 
maxima  in  the  amplitude.  However,  after  the  initial  complete  rotations,  the 
relative  phase  appears  locked  in  the  interval  it  <  4>^  -  ^  <  3n/2,  and 
accordingly,  the  frequency  of  the  Jl=l  mode  oscillates  about  that  of  the  £=2 
dominant  mode. 

Before  concluding,  let  us  point  out  that  in  all  two-mode  simulations,  those 
of  Table  2  and  those  not  presented  here,  only  one  mode  seems  to  dominate.  This 
mode  is  not  necessarily  the  one  corresponding  to  the  highest  amplitude  or 
efficiency.  For  example,  the  amplitude  and  efficiency  of  the  dominant  £= 2 
mode  in  cases  1  and  5  are  significantly  below  the  corresponding  values  of  the 
£=1  mode,  if  the  excitation  of  the  latter  is  assured  through  priming,  as  is 
done  is  cases  8  and  9.  The  opposite  holds  true  in  a  comparison  of  case  10  to 
case  7.  On  the  other  hand,  simulations  of  a  long  cavity,  in  which  many  modes 
are  resonant  with  the  beam,  ranging  from  £=1  to  $,=10,  have  shown  that  many  of 
these  modes  are  simultaneously  excited  to  comparable  amplitudes  and  efficiencies. 
Such  situations  have  not  been  adequately  analyzed  at  this  time,  and  accordingly 
we  choose  not  to  present  them  here. 


9.  DISCUSSION 


In  this  work  we  have  developed  the  self-consistent  formalism  for  the 
interaction  of  the  eigenmodes  with  the  electron  beam  in  an  overmoded  gyrotron 
cavity.  The  resulting  equations  have  been  applied  in  various  simplified 
situations,  both  analytically  and  numerically,  in  order  to  establish  their 
correctness,  to  extend  earlier  results,  to  interpret  the  effects  of  an  axial 
nonuniformity  and  to  initiate  the  investigation  of  the  multi-mode  interactions. 
The  approach  we  have  adopted,  discussed  in  Sec.  3,  consists  of  an  integration 
over  the  fast-time  scale,  which  includes  variations  at  the  wave  and  the 
electron  cyclotron  frequencies,  resulting  in  equations  depending  only  on  first 
order  parameters  describing  the  amplitude  of  the  cavity  fields  and  the  non¬ 
uniformities,  and  varying  only  at  the  slow-time  scale  characterized  by  the 
bounce  frequency  and  the  relative  frequencies.  The  derivation  of  these 
equations  was  discussed  in  Sec.  4.  In  addition  to  self-consistency,  Eqs.  (10) 
and  (12)  have  the  advantages  of  applicability  to  analytical  investigations  and 
of  high  efficiency  when  used  in  numerical  simulations. 

The  correctness  of  the  equations  was  established  in  Sec.  5,  and  the 
accuracy  of  their  numerical  implementation  was  demonstrated  in  Sec.  6,  by 
applying  a  simplified  version  of  the  equations  to  study  the  evolution  in  a 
rectangular  waveguide.  Our  analysis  has  enabled  us  to  derive  general  expres¬ 
sions  for  the  dispersion  relation  and  for  the  bounce  frequency  for  the  case 
of  arbitrary  variations  of  the  wave  in  both  transverse  directions.  These 
results  are  represented  by  Eq.  0 3)  (or  Eq.  (15)  for  conditions  well  above 
threshold)  and  Eq.  (18).  In  addition,  a  dimensionless  parameter.  A,  was 
obtained  In  Eq,  (17),  which  is  useful  to  determine  the  saturation  level  of  the 
instability,  and  simplified  equations  were  obtained  for  the  particle  dynamics 
and  wave  evolution  at  saturation  (see  Eqs.  (19)  and  (21)). 


The  subject  of  external  axial  nonuniformities  was  considered  in  Sec.  7. 

The  equations  of  motion  were  simplified  to  the  pendulum  equation,  Eq.  (22), 
with  an  external  torque,  given  by  R  in  Eq.  (23)  (or  Eqs.  (26)  or  (27)) 
representing  the  effects  of  the  nonuni fonni ties.  Based  on  the  corresponding 
electron  trajectories,  the  efficiency  enhancement  was  interpreted  when  R  >  0, 
while  the  possibility  of  similar  effects  was  discussed  for  the  case  R  <  0, 
under  the  appropriate  conditions.  It  was  also  shown  that  a  positive  efficiency 
should  be  expected  also  in  the  case  of  a  beam  with  substantial  thermal  spread, 
both  when  R  <  0  and  R  >  0,  under  conditions  roughly  representative  of  an 
oscillator  and  amplifier,  respectively.  Finally,  the  equivalence  of  a  multi¬ 
tude  of  nonuniformities  was  shown.  Such  nonuniformities  include  a  gradient 
of  the  external  magnetic  field,  a  gradient  in  the  cross  section,  the  presence 
of  an  axial  static  electric  field  and  the  introduction  of  a  dielectric  material 
with  a  dielectric  constant  which  depends  on  position. 

Some  introductory  investigation  of  two-mode  competition  has  been  presented 
in  Sec.  8.  It  was  seen  that  one  mode  eventually  dominates  the  spectrum  and 
effectively  suppresses  the  excitation  of  the  second  mode,  regardless  of  whether 
or  not  the  unamplified  mode  would  have  reached  higher  power  levels,  were  it 
permitted  to  grow  alone.  This  behavior  is  similar  to  that  observed  by 
Dialetis  and  Chu  (21).  Analytical  expressions  for  the  saturation  amplitude 
and  efficiency  were  obtained,  see  Eqs.  (28)  and  (29),  and  satisfactorily 
verified  by  the  simulations.  Finally,  it  was  shown  that  the  asymptotic 
behavior  of  the  system  is  characterized  by  amplitude  and  frequency  oscillations 
at  the  beat  frequency.  It  can  easily  be  seen  that  such  oscillations  correspond 
to  a  weak  excitation  of  a  discrete  sideband  spectrum  with  separation  equal  to 
the  frequency  difference  of  the  two  modes.  It  is  possible  that  the  simultaneous 
excitation  of  many  modes  in  a  highly  overmoded  cavity  is  related  to  the 
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corresponding  denseness  of  this  spectrum  parts  of  which  might  well  resonate 
effectively  to  different  eigenmodes. 

At  this  point  it  is  appropriate  to  present  a  discussion  on  the  limitations 
of  the  test  particle  analysis  by  Dialetis  and  Chu  (21).  One  limitation  in 
their  approach  is  related  to  the  absence  of  self-consistency,  caused  by  the 
omission  of  the  reactive  component  of  the  current.  This  suppresses  any  shift, 
A,  of  the  frequency  from  its  empty  cavity  value,  causing  incorrect  conclusions 
to  be  drawn  when  this  shift  is  important.  Such  a  situation  would  arise  when 
the  frequency  shift  is  comparable  to,  or  larger  than,  either  the  frequency 
mismatch,  or  the  inverse  transit  time,  that  is,  either  A/e  £  1  or  :  A  >  1 , 
since  in  such  cases  one  would  obtain  inaccurate  evaluations  of  the  motion  of 
the  trapped  or  untrapped  particles,  respectively.  In  general,  however, 

A~ui/Q,  as  can  be  seen  from  Eq.  (12),  hence  the  parameters  they  consider  satisfy 
A/e  -  tA—0.1  ,  and  no  danger  of  substantial  error  appears  present.  More 
serious  could  be  their  assumption  on  the  existence  of  a  steady  state.  As  has 
been  shown  in  our  simulations,  the  equilibrium  state  is  not  steady,  but 
oscillatory.  These  oscillations  may  be  very  important  when  more  than  two 
modes  are  simultaneously  present,  because  of  the  possibility  of  synchronism 
between  parts  of  the  discrete  spectra,  which  such  oscillations  represent. 

In  conclusion,  let  us  discuss  various  extensions  and  modifications  to  our 
work.  One  desirable  modification  appears  to  be  that  associated  with  the 
adoption  of  a  cylindrical  geometry,  which  is  more  compatible  with  the  experi¬ 
mental  arrangements.  No  conceptual  difficulty  is  expected,  with  the  method 
of  Sec.  3,  including  Eq.  (6)  being  directly  applicable.  The  eigenmodes  will 
have  a  different  transverse  dependence,  and  the  only  added  feature  will  be  the 
application  of  Graf's  addition  theorem  In  the  evaluation  of  the  fields  along 
the  zero  order  trajectories,  e.g.. 


Jm(knr)  cos  (me)  =  £  J,n+q(kn^  cos  [m^e  "  8)  +  q(ir  +  ^  -  §)]  , 


(30) 


where  tan8  =  y/x  and  tane  =  y/x,  while  here  n  and  m  are  the  radial  and 
azimuthal  harmonic  numbers. 


A  second  possible  extension,  independent  of  the  choice  of  rectangular  or 
cylindrical  geometry,  involves  the  introduction  of  alternate  axial  nonuniformi- 
ties.  The  additional  force  due  to  an  axial  electric  field  is  a  trivial  exten- 
tion,  while  the  inertial  forces  produced  by  a  variation  of  the  cross  section 
and  the  introduction  of  a  dielectric  sleeve  with  position  dependent  dielectric 
constant  can  be  accounted  for  by  assigning  a  WKB-type  axial  field  dependence, 
of  the  form  sin  k^z,  where  k^  is  now  a  function  of  z,  calculated  in  terms  of 
the  frequency  and  the  local  value  of  the  transverse  wave  number.  In  addition, 
and  regardless  of  the  above,  calculations  have  to  be  performed  to  extend  the 
definition  of  R  (Eq.  (23))  to  the  high  efficiency  regime,  to  investigate  the 
dependence  of  efficiency  enhancement  on  R,  and  to  assess  the  effects  of  the 
direction  of  the  nonuniformity  and  the  possibility  of  obtaining  acceptable 
values  for  the  efficiency  in  spite  of  a  thermal  spread,  as  discussed  in  Sec.  7. 

Finally,  the  introductory  results  of  Sec.  8  are  expected  to  be  useful  in 
various  aspects  of  the  study  of  multi -mode  competition.  An  interesting  first 
step  appears  to  be  the  investigation  of  particle  trajectories  in  a  two-mode 
field,  with  a  steady  large  amplitude  dominant  mode  accompanied  by  an  oscilla¬ 
tory  low  amplitude  secondary  mode,  to  establish  criteria  for  the  stability 
of  such  a  configuration.  The  following  step  would  be  to  study  the  effects  of 
coupling  between  harmonics  of  the  beat  frequencies  in  a  highly  overmoded 
cavity,  In  which  more  than  two  resonant  modes  are  present.  In  addition,  the 
competition  of  modes  with  different  tranverse  structure  should  be  considered. 
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Clearly,  all  these  subjects  are  to  be  pursued  primarily  numerically,  with 
supporting  analytical  considerations to  the  extent  it  is  feasible. 
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